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Abstract 


This study considers whether spurious peaks in nitrate ions in snow sampled at Summit, 
Greenland from August 2000 to August 2002 are related to solar proton events. After 
identifying tropospheric sources of nitrate on the basis of correlations with sulfate, 
ammonium, sodium, and calcium, we use the three-dimensional global Whole 
Atmosphere Community Climate Model (WACCM) to examine unaccounted for nitrate 
spikes. Model calculations confirm that solar proton events significantly impact HO x , 
NO x , and 0 3 levels in the mesosphere and stratosphere during the weeks and months 
following the major 9 November 2000 solar proton event. However, SPE-enhanced NO y 
calculated within the atmospheric column is too small to account for the observed nitrate 
ion peaks in surface snow. Instead, our WACCM results suggest that nitrate spikes not 
readily accounted for by measurement correlations are likely of anthropogenic origin. 
These results, consistent with other recent studies, imply that nitrate spikes in ice cores 
are not suitable proxies for individual SPEs and motivate the need to identify alternative 


proxies. 
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1. Introduction 

20 Identifying the impact of solar particle storms on the atmosphere remains fundamental 
in understanding the Sun's influence on Earth's climate [ Gray et al., 2010; National 
Research Council, 2012]. High-energy particles from these solar events increase odd 
nitrogen and odd hydrogen, catalytically destroying ozone and thereby potentially 
impacting climate through the chemistry, radiative budget, and dynamics of the upper 
25 atmosphere [e.g., Randall et al., 2005; Jackman et al., 2008]. In addition, these space 
weather events have the potential to disrupt power grids, communications technology, 
and spacecraft [National Research Council, 2008]. 

Direct observations of solar energetic particle events have only been available since 
the mid-20th century. A broader understanding of the potential frequency and intensity of 
30 these events requires a more extensive record of historical occurrences, motivating the 
search for indirect proxy evidence [Schrijver et al., 2012]. There is extensive analytical 
and predictive research using nitrate variability in polar ice cores as a proxy for solar 
energetic particle events [e.g., Zeller and Parker, 1981; Dreschhoff and Zeller, 

1990; McCracken et al., 2001a; Shea et al, 2006; Kepko et al., 2009]. However, this 
35 relationship has been questioned, particularly with regard to the short timescales 

associated with individual events [e.g., Legrand and Delmas, 1986; Wolff et al., 2008; 
Wolff et al., 2012]. Contemporary progress toward predicting space weather urgently 
awaits the resolution of whether or not nitrate ion spikes in ice cores can be used to infer 
past events [e.g., Barnard et al., 20 1 I ; Riley, 2012]. 

40 Zeller and Parker [1981] associated nitrate levels with solar activity through the 

correlation of nitrate ions (NO3") in Antarctic ice cores with cosmogenic carbon isotopes 
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( 14 C) in tree rings. Statistical correlation studies confirm the covariance between nitrate 
ions and cosmogenic radionuclides 14 C and 10 Be on centennial to millennial timescales 
[e.g., McCracken et al., 2001b; Traversi et al., 2012; Ogurtsov and Oinonen, 2014]. 

45 However, attempts to find a correlation between nitrate ions and solar variability on 
smaller timescales, such as the 1 1-year solar cycle or individual events, have been 
unsuccessful [e.g., Legrand and Delmas, 1986; Legrand and Kirchner, 1990; Legrand et 
al, 1996; Traversi et al., 2012]. Instead, the variability of nitrate ions in polar ice is 
attributable to lightning from lower latitudes and downward transport from the lower 
50 stratosphere [Legrand and Delmas, 1986; Legrand et al., 1989; Legrand et al., 1996], 
with potential contributions in the Arctic from anthropogenic pollution [Mayewski et al., 
1990] and biomass burning [e.g., Whitlow et al., 1994; Dibb and Jafrezzo, 1997; 

Savarino and Legrand, 1998]. 

Efforts to attribute sharp nitrate ion peaks in ice cores to individual solar proton events 
55 are exemplified by the Zeller and Dreschhoff [1995] analysis of 8000 years of nitrate data 
from the GISP2-H Greenland ice core and the estimated cumulative probabilities of solar 
event occurrences by McCracken et al. [2001a]. Palmer et al. [2001] provide 
observational evidence of a small statistical background increase in nitrate ions as a result 
of individual SPE events. However, theoretical considerations and model simulations 
60 suggest that enhancements of nitrogen species from individual solar events are unlikely to 
produce sharp peaks in nitrate ions at the surface, given the slow rate of vertical transport 
in the stratosphere, horizontal mixing and dilution to lower latitudes, and diabatic 
recirculation [e.g., Legrand et al., 1989; Legrand and Kirchner, 1990]. In addition, 

Weller et al. [2011] find no indication of individual solar events in 25 years of 
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65 atmospheric aerosol measurements in Antarctica. Wolff et al. [2008; 2012] are able 

instead to attribute nitrate ion peaks to tropospheric sources using correlations among a 
suite of ions in surface snow. 

This paper presents a case study adding to a growing literature that challenges the 
validity of using nitrate spikes in ice as proxies for individual solar energetic particle 
70 events. A three-dimensional global simulation of the 9 November 2000 solar proton event 
using the Whole Atmosphere Community Climate Model (WACCM) combined with 
daily samples of nitrate ions in surface snow at Summit, Greenland demonstrates the 
difficulty in generating impulsive enhancements in nitrate deposition in snow and ice 
from a single event. The results of this study, however, reiterate the significant influence 
75 of solar energetic particles on the chemical composition of the Arctic polar stratosphere 
and mesosphere, especially levels of odd nitrogen and ozone, encouraging the search for 
robust proxies other than nitrate to determine the frequency and intensity of historical 
solar events. 

so 2. Methods 

2.1. Schematic of Solar Protons Events and Nitrate Precursors 

Solar energetic particles precipitating into our atmosphere include electrons, protons, 
and more massive ions. Our study focuses on solar proton events (SPEs), as solar protons 
uniquely possess sufficient energies to penetrate and modify the chemistry of the 
85 mesosphere and stratosphere, potentially leading to chemical signatures within the 

troposphere and at the surface. Figure 1 presents a schematic of the perturbations caused 
by solar energetic protons within the polar atmosphere. High-energy protons from solar 
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flares and coronal mass ejections precipitate over the polar caps, commonly extending to 
geomagnetic latitudes greater than 60 degrees [Smart and Shea, 1994]. These high- 
90 energy protons, along with collisional secondary electrons, ionize and dissociate 

molecular nitrogen and oxygen, resulting in the formation of odd hydrogen (HO x = H + 
OH + HO 2 ) and reactive odd nitrogen (NO x = N + NO + NO 2 ) [e.g., Crutzen et al., 1975; 
Jackman et al., 1980; Solomon et al., 1981]. 

Enhancements of HO x cause short-lived catalytic ozone destruction in the mesosphere 
95 and upper stratosphere during an SPE and for a few days following the event [e.g., 
Solomon et al., 1981; Jackman et al., 2008; Damiani et al., 2010]. NO x has a short 
lifetime in the upper mesosphere and thermosphere, but the lifetime increases to months 
during polar night in the stratosphere. During winter, when downward transport within 
the isolated polar vortex is strong and photochemistry is limited, NO x produced by SPEs 
100 in the mesosphere and upper stratosphere can be transported to the middle and lower 
stratosphere [Lopez-Puertes et al., 2005; Randall et al., 2005; Jackman et al., 2009; 
Randall et al., 2009]. Loss of O 3 from the oxidation of SPE-enhanced NO x mainly results 
in nitric acid (HNO 3 ) and dinitrogen pentoxide (N 2 O 5 ). These species eventually mix 
with the large pool of total odd nitrogen (NO y = N + NCE + NO 3 + 2N 2 0s + HNO 3 + 

105 HO 2 NO 2 + (HONO) + CIONO 2 + BrONCE) in the lower stratosphere, a background 
reservoir produced primarily by the oxidation of nitrous oxides (N 2 0) emitted at the 
surface [Vitt and Jackman, 1996]. With the return of sunlight to the polar region, NO x can 
more effectively destroy stratospheric ozone through catalytic reactions. 

Total odd nitrogen species such as HNO 3 , H0 2 N0 2 , C10N0 2 , and N 2 0 5 serve as 
110 precursors to soluble nitrate ions (NO 3 ") deposited in snow. Nitrate ions from HNO 3 , and 
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to a less extent HO2NO2, can reach the surface through stratosphere-troposphere 
exchange followed by wet deposition. In the winter polar vortex, HNO3 also condenses to 
form polar stratospheric cloud (PSC) particles, which can gravitationally sediment to the 
troposphere if the particles grow large enough. In addition, heterogeneous reactions of 
115 N2O5 + H2O and CIONO2 + HC 1 occur on the surfaces of PSCs to form HNO3, 

condensing and sedimenting out of the stratosphere. On the basis of inferences from 
nitrate layers in ice core data, McCracken et al. [2001a] and Shea et al. [2006] suggest 
that the deposition of nitrate to the surface resulting from SPEs occurs 2-6 weeks 
following each event. This contradicts calculations by Legrand et al. [1989] that suggest 
120 a two-year transport time of enhanced odd nitrogen from the upper stratosphere to lower 
stratosphere, a consequence of horizontal diffusion to lower latitudes and vertical diabatic 
recirculation. 

A viable mechanism for SPE nitrate precursors to progress from the upper atmosphere 
to the surface snow within a 2-6 week timeframe would require: 1) rapid downward 
125 transport from the upper stratosphere; 2) levels of NO y produced by SPEs high enough to 
compete with the background reservoir of NO y in the lower stratosphere; and 3) a 
mechanism for quickly depositing nitrate from the lower stratosphere to the surface. High 
temporal resolution would be necessary to identify these events at the surface, requiring 
ice core sampling techniques such as continuous flow analysis (CFA) [e.g., Sigg et al., 

130 1994; Roethlisberger et al., 2000; Kepko et al., 2009] or high frequency sampling of 

surface snow or ambient air [e.g., Wolff et al, 2008; Weller et al., 2011]. 

Using daily measurements of surface snow at Summit, Greenland, the investigation 
described in this paper considers much higher resolution nitrate variability than prior ice 
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cores studies, resolving nitrate deposition on timescales that should be able to capture an 
135 individual SPE. Although there are larger SPE events during meteorological periods of 
stronger downward transport, we target the SPE of 9 November 2000 because daily 
surface snow measurements following this event are complete enough to infer alternative 
tropospheric sources of nitrate ion spikes through correlation analysis with sulfate, 
ammonium, sodium, and calcium as discussed in Section 3.1. The 9 November 2000 
140 event is the sixth largest SPE in the last 50 years with respect to the calculated production 
of NO y , four times weaker than the largest event on 19-27 October 1989 [. Jackman et al., 
2008]. 

2.2. The Whole Atmosphere Community Climate Model 

145 The Whole Atmosphere Community Climate Model (WACCM) is a component of the 
Community Earth System Model (CESM) at the National Center for Atmospheric 
Research (NCAR). Model documentation is available at the CESM website 
(www2.cesm.ucar.edu) . Simulations in this study use version cesml.0.5 (WACCM4), 
with active atmospheric and land models, prescribed ice, and fixed ocean (specified sea 
150 surface temperatures). The WACCM atmospheric component of the model combines the 
Community Atmosphere Model (CAM5), the Thermosphere-Ionosphere-Mesosphere- 
Electrodynamics General Circulation Model (TIME-GCM), and the Model for OZone 
and Related chemical Tracers (MOZART) to simulate dynamics and chemistry from the 
surface to the lower thermosphere [Garcia et al, 2007; Kinnison et al., 2007; Emmons et 
155 al., 2010; Marsh et al., 2013; Neale et al., 2013]. 
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The chemical solver and reaction rates are based on MOZART chemistry [ Kinnison et 
al., 2007; Emmons et a/., 2010]. Observed solar spectral irradiance and geomagnetic 
activity force the heating and photolysis rates [Marsh et al., 2007]. The chemical 
mechanism includes 59 species and involves reactions of O x , NO x , HO x , C10 x , and BrO x 
160 chemical families as well as methane and carbon monoxide oxidation. The mechanism 
contains heterogeneous reactions on stratospheric aerosols, including liquid sulfate 
aerosols along with nitric acid trihydrate (NAT), supercooled ternary solution (STS), and 
water ice associated with PSCs [Kinnison et al, 2007]. Concentrations of longer-lived 
greenhouse gases and halogen species are specified from observations [Garcia et al., 

165 2007]. Surface emissions are represented by flux boundary conditions associated with the 

most recently available compilation supported by WACCM and described by Lamarque 
et al. [2012]: anthropogenic emissions from POET (Precursors of Ozone and their Effects 
on the Troposphere) [Granier et al, 2005] and REAS (Regional Emissions inventory in 
ASia); monthly biomass burning emissions from GFED-v2 (Global Fire Emissions 
no Database) [van der Werf et al., 2006]; and biogenic, soil, ocean, and volcanic emissions 
from POET and GEIA (the Global Emissions Inventory Activity). These inventories 
correspond to emissions used in the Coupled Model Intercomparison Project phase 5 
(CMIP5) [Taylor et al., 2012]. 

The WACCM simulations used here have a resolution of 1.9° latitude, 2.5° longitude, 
175 and 88 vertical layers extending from the surface to approximately 140 km. The model 
chemistry applies 30-minute time steps. We use the specified dynamics version of 
WACCM (SD-WACCM) in which meteorology is forced by NASA's Modern Era 
Retrospective analysis for Research Applications (MERRA) fields [Rienecker et al., 
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2011]. WACCM achieves this forcing by relaxing horizontal winds and temperatures to 
180 MERRA fields from 0 to 40 km. The model is free running above 50 km, with a linear 
reduction of forcing between 40 and 50 km. 

WACCM modelers have participated in a series of intercomparisons and validation 
studies, most recently CMIP5 [Taylor et al., 2012]. Funke et al. [2011] include WACCM 
in an intercomparison focused on the SPEs of October 2003, validating results with 
185 observations from the Michelson Interferometer for Passive Atmospheric Sounding 

(MIPAS) on Envisat. The authors show a general agreement between model calculations 
and measurements, including O 3 loss within 5% and model NO y enhancement within 
30% at 1 hPa. 

Previous WACCM simulations studying SPEs include Jackman et al. [2008; 2009; 

190 2011] and Funke et al. [2011]. The WACCM simulations presented in this paper are 

motivated in part by Figure 6 from Jackman et al. [2009], showing a tongue of enhanced 
NO y and decreased O3 extending into the northern polar lower stratosphere following the 
9 November 2000 SPE. The present simulations differ by: 1) forcing the model with 
MERRA re-analysis meteorological fields specific to the 2000-2001 time period; 2) using 
195 a more recent version of WACCM with improvements including parameterizations for 
wave interactions resulting from turbulent mountain stress; and 3) using a higher 
horizontal resolution to better resolve transport (1.9° latitude by 2.5° longitude in contrast 
to the Jackman et al. [2009] 4° latitude by 5° longitude). 


200 2.3. Modeling the November 9, 2000 Solar Proton Event 
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Following a coronal mass ejection on 8 November 2000 (23:26 UTC), the GOES8 
satellite measured a solar proton event (SPE) beginning on 8 November (23:50 UTC) and 
reaching a maximum on 9 November (15:55 UTC). A smaller event occurred in late 
November, peaking on 26 November (20:30 UTC). The proton flux for energies greater 
205 than 10 MeV reach a maximum of 14,800 pfu (pfu = proton flux units = 

particles-sr ' -cm' -s' ) for the 9 November event and 940 pfu for the 26 November event 
( http://www.swpc.noaa.gov) . 

The WACCM model assumes a unifonn distribution of proton flux at geomagnetic 
latitudes greater than 60 degrees, using calculations of daily-averaged ion pair production 
210 rates as a function of pressure based on GOES 8 proton flux measurements [. Jackman et 

al., 1980; 2005; 2008]. WACCM includes HO x production rates as a function of altitude 
and ion pairs using a table from Jackman et al. [2005] based on the dissociation of O 2 
followed by water cluster ion formation and neutralization described by Solomon et al. 
[1981]. NO x formation results from the ionization and dissociation of N 2 and 0 2 , 

215 producing ~1 .25 N per ion pair [ Porter et al., 1976] partitioned into 45% ground state 
N( 4 S) and 55% excited state N( 2 D). The excited state N( 2 D) determines net NO x 
production [Rusch et al, 1981]. Tabulated calculations and detailed descriptions of NO x 
and HO x production rates by SPEs are available at the SPARC/SOLARIS website 
(http://sparcsolaris.geomar.de/input_data.php). 

220 WACCM calculations of gas phase chemical species include loss rates for dry 

deposition, wet deposition, and heterogeneous reactions on stratospheric aerosols. Dry 
deposition follows a resistance approach, dependent on land cover type and surface 
roughness. WACCM calculates the dry deposition flux for a given species as the product 
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of deposition velocities (varying by time and horizontal grid location) and concentrations 
225 at a reference height (~10 m) above the surface. 

There are two wet deposition schemes available for WACCM: the traditional 
MOZART scheme [Rasch et al, 1997; Horowitz et al., 2003] and a more extensive 
scheme detailed by Neu and Prather [2012]. The MOZART scheme addresses convective 
updrafts, in-cloud nucleation scavenging (rainout), below-cloud impaction scavenging 
230 (washout), and evaporation in clear ambient air. Precipitation rates, cloud water content, 
and cloud fractions are taken from meteorological data fields, and effective Henry's law 
coefficients are prescribed. The Neu and Prather [2011] wet deposition scheme includes 
a more resolved sub-grid scale treatment of cloud overlap and a burial method for the 
uptake of soluble gases on ice, resulting in slower uptake of HNO3 by ice and snow. This 
235 study uses the MOZART scheme in an effort to provide an upper limit for wet deposition 
of HNO3, consistent with the search for the maximum potential deposition of nitrate 
following SPEs. 

WACCM treats wet deposition as a first order gas-phase loss process at the end of 
each time step. The model does not explicitly account for the accumulation of condensed 
240 species within cloud droplets or aqueous chemistry, nor does it archive soluble ions 
deposited to the surface. The results presented in this paper estimate the variability in 
nitrate deposition by integrating gas-phase loss through wet deposition over the total 
atmospheric column and then dividing by precipitation amounts, filtering model output 
when precipitation values are extremely low to prevent anomalous spikes. These 
245 calculations are expected to overestimate the magnitude of nitrate ions deposited to the 
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snow but allow a comparison of modeled variability in nitrate deposition with observed 
variability in surface snow measurements. 

WACCM also calculates the uptake of HNO3 by stratospheric aerosols and the 
gravitational settling of stratospheric cloud particles. Combining the sedimentation of 
250 condensed phase HNO3 with the heterogeneous conversion of N2O5 to HNO3 on polar 
stratospheric clouds allows WACCM to calculate the contribution of nitrate to the 
troposphere through irreversible denitrification of the stratosphere. 

2.4. Observations of Nitrate in Surface Snow at Summit, Greenland 

255 Summit Station, originally the site of the Greenland Ice Sheet Project 2 (GISP2), is 
located in the middle of the Greenland ice sheet, approximately 3200 meters above sea 
level at 72° 34' N latitude and 38° 29' W longitude. Atmospheric measurements have been 
made at the site since 1989, with a variety of intensive measurements throughout the 
1980s and 1990s, culminating in a full suite of continuous measurements beginning in 
260 2 0 03 as part of the Arctic Observing Network (http://www.geosummit.org) . 

Year-round daily measurements of soluble ion content in surface snow are available 
during 1997-1998 and from August 2000 to August 2002 [Dibb et al., 2007]. Vertical 
profiles from monthly one-meter snow pits accompany these measurements. Sampled 
ions include sodium (Na + ), ammonium (NH 4 + ), potassium (K + ), magnesium (Mg 2+ ), 

265 calcium (Ca 2+ ), chloride (Cf), nitrate (N0 3 ‘), and sulfate (S0 4 2 '). Sampling procedures, 
transport, and ion chromatography techniques and uncertainties are described in detail by 
Dibb et al. [2007]. Briefly, snow was sampled daily from an areas upwind of Summit 
camp to avoid contamination from local sources. Samples from the uppermost 
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stratigraphic layer were collected in groups of three adjacent replicates, with companion 
270 blanks generated every nine samples. Samples remained at -20 C until analysis, melted in 
small batches and immediately analyzed by ion chromatography. 

The use of chemical profiles to infer historical atmospheric conditions such as the 
influence of solar proton events depends on preservation within the snowpack. The 
analysis of surface snow and snow pit data reveals that levels of most ions are 
275 preserved (at approximately 80%) in layers within the snow pits [Dibb et al., 2007]. The 
exception is nitrate, where post-depositional processes may significantly modify nitrate 
ion concentrations [Dibb et al., 2007]. Evidence of post-depositional modifications have 
been observed at both poles [e.g., Dibb and Whitlow, 1996; Dibb and Jaffrezo, 1997; 
Rothlisberger et al., 2002], with significant differences expected among sites as a result 
280 of the physical and chemical environments. These modifications involve snow-air 
exchange attributed to photolysis, sublimation and condensation associated with the 
growth of snow grains, and the uptake and release of volatile species such as nitric acid 
[e.g. Legrand et al., 1 996; Sturm and Benson, 1997; Honrath et al., 1999; Dibb et al., 
2002; Gr annas et al., 2007]. Traversi et al. [2012] suggest that ideal conditions for 
285 preserving nitrate in firn with respect to post-depositional effects involve accumulation 
rates above 50 mm water equivalent per year. While sites in Greenland such as Summit 
are well above this threshold, many Antarctic plateau sites fall below the limit, including 
Vostok, and Dome C [ Traversi et al., 2012 and references therein]. In addition, long 
periods between snowfalls may allow significant post-depositional processing near the 


290 surface. 
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Post-depositional modifications also depend on the chemical composition of air near 
the surface, the composition of the snow (particularly with regard to acidity), and the 
altitude and temperature of the site [e.g., Rothlisberger et al., 2000; Burkhart et ah, 

2004]. Recent measurements of oxygen isotopes of nitrate at Summit suggest that post- 
295 depositional photolysis of nitrate in surface snow is limited (~2% in summer) [Fibiger et 
al., 2013]. Furthermore, in the absence of solar radiation, we would not expect post- 
depositional loss process involving photochemistry and thermal desorption to extensively 
reduce concentrations during the polar winter. However, modifications of nitrate ion 
levels may result from the migration, diffusion, and chemistry associated with the 
300 physical evolution of the snowpack, particularly with regard to the metamorphism and 
compaction of grains [ Bartels-Rausch et al., 2012]. 

The seasonal cycles of long-lived chemical tracers, greenhouse gases, and non- 
methane hydrocarbons demonstrate that this site is representative of well-mixed, remote 
Arctic air, with most organic compounds in their final stage of oxidation [ Dibb et ah, 

305 2007]. Kahl et al. [1997] show westerly flow dominating at Summit, with wintertime 

trajectories at 500 hPa identifying long-range transport from Asia or Europe and 
trajectories at 700 hPa suggesting less rapid flow from North America. 

Summit receives ~65 cm of snow (—24 cm water equivalence) per year, with snow 
accumulation reaching a maximum in spring and a minimum in winter 
3io [Dibb and Fahnestock, 2004]. Periods of fresh snow are relatively infrequent during the 
2000-2001 winter: six events in November 2000 (0.72 cm water equivalence); two events 
in December 2000 (0.29 cm water equivalence); and two events in January 2001 (0.56 cm 
water equivalence). In addition to fresh snow events, field notes indicate frequent periods 
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of thick fog and rime, ice-crystals and diamond dust as well as high winds and blowing 
315 snow strong enough to obscure visibility. 

3. Results and Discussion 

3.1. Correlations Between Nitrate Ions and Other Ions in Surface Snow at Summit 

Concurrent measurements of nitrate (NO3') with other ions such as sulfate (SO4 2 '), 

320 calcium (Ca 2+ ), sodium (Na + ), and ammonium (NH 4 + ) provide a means for attributing 

nitrate spikes in snow to tropospheric sources, either by direct association with industrial 
pollution and biomass burning plumes from Europe and North America or from increased 
nitrate deposition involving particulates such as sea salt and dust. Mavewski et al. [1990] 
credit a rising trend in both NO3' and SC>4 2 ‘ in Greenland to industrial continental 
325 pollution. Anthropogenic pollution not only is characterized by high levels of HNO3 but 
also is the dominant source of sulfate in comparison to sea salt, dust, volcanoes, and 
biogenic emissions. Chemical signatures of continental biomass burning include elevated 
levels of NH4 + and formate in Greenland ice, with concurrent enhancements of NO3' 
present in some events and absent in others [e.g., Legrand et al., 1992; Whitlow et al., 

330 1994; Legrand and de Angelis, 1996; and Savarino and Legrand, 1998]. Dibb et al. 

[1996] and Dibb and Jafrezzo [1997] confirm biomass burning in Northern Canada as a 
source of enhanced NH 4 + , carboxylic acids, and nitrate at Summit using back- trajectories 
in conjunction with with atmospheric and snow measurements. 

Observations of atmospheric aerosols at Summit near ground level (1.5 m) show 
335 minimal concentrations of particulate nitrate, often below the detection limit of 

instruments, suggesting that the majority of atmospheric nitrate in snow originates from 
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gas phase HNO 3 [Dibb et al., 1994; Dibb et al., 1998]. Contemporaneous measurements 
of gaseous HNO 3 , atmospheric aerosol nitrate, and nitrate ions in snow during the 

■j 

summer of 1993 show mean HNO3 concentrations of 0.9±0.6 nmol/m an order of 
340 magnitude larger than mean aerosol nitrate concentrations 0.06±0.6 nmol/m [Dibb et al, 
1994]. N 2 O 5 could potentially be an additional source of nitrate ions in snow, particularly 
in winter, when N 2 O 5 accumulates in the absence of photochemistry and subsequently 
reacts heterogeneously with the snow [Huff et al., 2011]. 

Wolff et al. [2008] attribute correlations between NO3' and Na + at Halley, Antarctica 
345 to increased rates of conversion of gaseous to aerosol nitrate on coarse sea salt aerosols or 
salty snow surfaces. The authors propose the potential for similar correlations involving 
NH 4 H and Ca 2+ at inland polar sites like Summit. Dibb et al. [2007] identify a spring 
maximum in crustal dust (Ca ) at Summit, most likely long-range transport originating 
from Asia [Dibb et al., 2003] as well as a summer peak in NH [Dibb et al., 2007], 

350 characteristic of biomass burning from North America and Europe. Wolff et al. [2008] 
recommend that future analyses of nitrate ions in snow and ice use ion correlations to 
screen for known tropospheric sources before searching the residual nitrate data for 
potential SPEs. This study adopts such a technique, identifying tropospheric sources 
through ion correlations at Summit from 2000-2002 followed by WACCM model 
355 simulations to investigate nitrate spikes not readily attributable to tropospheric sources. 

Figure 2 presents time series plots of NO 3 ' paired with NH 4 + , Na + , Ca 2+ , and S0 4 2 ' at 
Summit along with solar proton fluxes from GOES 8 . Visual examination shows that 
nitrate spikes on 21 August 2000, 5 September 2000, 16 June 2001, and 4 July 2002 
correlate with NH 4 + , a tracer of biomass burning. High levels of the anthropogenic 
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360 pollutant tracer SOW accompany nitrate spikes on 25-30 November 2000, 26 June 2001, 
and 19 February 2002. Levels of Na + associated with sea salt are enhanced along with 
NO 3 ' on 12 January 2002 and 19 February 2002. This leaves four nitrate ion spikes that 
are not readily attributable to tropospheric sources: 22-24 November 2000, 13 December 
2000, 25 January 2001, and 22 October 2001. The GOES 8 proton fluxes indicate SPEs 
365 on 9 Nov 2000, 24 Sep 200 1 , 4 Nov 200 1 , and 22 Nov 200 1 . This study uses the 

WACOM model to analyze the three candidate nitrate spikes that occur during polar 
winter (22-24 November 2000, 13 December 2000, 25 January 2001), specifically 
searching for a potential mechanism to explain how a nitrate signal from the 9 November 
2000 SPE might impulsively arrive at the surface through chemistry, transport, and 
370 deposition processes within the winter polar vortex. 

3.2. Modeling the Impact of the November 9 SPE on HO x , NO x , and O 3 

Comparing WACCM simulations with and without SPEs isolates the impact of solar 
protons on upper atmospheric processes and identifies potential scenarios leading to 
375 enhanced nitrate deposition. The model is initially run from January through October to 
allow chemical species to reach equilibrium throughout the atmosphere, with solar 
protons (e.g., the 14 July 2000 "Bastille Day" solar proton event) included during this 
“spin-up” period. At the start of November, the solar proton flux is allowed to continue in 
a “with SPEs” simulation but is set to zero in a “no SPEs” simulation, thereby more 
380 effectively highlighting perturbations specific to the 9 November SPE. 

Figure 3 shows the time evolution of the vertical structure of modeled HO x , NO x , and 
O3 (in mole ratios) during November. Top plots show WACCM simulations with SPEs. 


17 

Middle plots show WACCM simulations without SPEs. Bottom plots show the SPE 
enhancements of HO x and NO x and the SPE reduction of O 3 . These contour plots 
385 represent average profiles calculated within the meandering polar vortex, the region most 
likely associated with strong polar winter descent. Potential temperature (0) serves as the 
vertical coordinate on the left axis, while the right axis indicates the approximate 
geopotential height (Z). The boundary of the polar vortex is objectively determined by 
identifying grid points within the stratosphere where scaled potential vorticity (sPV) 

390 calculated on isentropic surfaces exceeds 1.4xl0' 4 s' 1 [Dunkerton and Delisi, 1986; 

Brakebusch et a/., 2013]. Scaled potential vorticity retains the conservation properties of 
Ertel's potential vorticity on isentropic surfaces. However, by normalizing with respect to 
U.S. Standard Atmosphere [ COESA , 1976], sPV removes vertical disparities in potential 
vorticity caused by the exponential increase of potential temperature with decreasing 
395 pressure. The vortex is assumed to remain in a constant location above 2500 K (~55 km), 
as the sPV method no longer adequately delineates the vortex edge as a result of the 
temperature profile above the stratopause. 

Figure 3 highlights the significant impact of the 9 November SPE event on the 
chemistry of the Arctic mesosphere and stratosphere. Short-lived enhancements up to 100 
400 ppbv HO x occur above 80 km in the mesosphere, with enhancements of 0.1 to 1 

ppbv extending throughout the stratosphere during the days following the 9 November 
SPE (Figure 3a). SPE production of NO x peaks around 50 ppbv in the mesosphere, 
exceeding 30 ppbv throughout the upper stratosphere during the days following the event 
(Figure 3b). The NO x enhancements slowly descend within the stratosphere throughout 
405 the month, maintaining levels an order of magnitude above background. Losses of O 3 up 
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to 800 ppbv occur above 70 km immediately following the 9 November SPE (Figure 3c). 
A 30 to 40% depletion of O 3 (> 100 ppbv) occurs throughout the stratosphere, most likely 
associated with the short-lived enhancement of HO x . O3 reduction continues in the 
stratosphere throughout the rest of November, consistent with the descent of SPE- 
410 enhanced NO x . The smaller SPE at the end of November is also evident in the plots for 
all three species. 

Figure 4 presents the vortex-averaged SPE enhancement of NO x and reduction of O 3 
from November through March. NO x enhancements of 10 to 40 ppbv (10 to 20 times 
background levels) descend from the upper stratosphere to the middle stratosphere 
415 throughout November and December at a rate of ~10 km/month. By January, increases of 
3 to 5 ppbv (2 to 3 times background levels) persist from 30 to 35 km, with remnant 
enhancements continuing into spring. 

Figure 4 supports the supposition that the enhancement of NO x during the months 
following the 9 November SPE drives longer-lived destruction of O3 in the stratosphere. 
420 The reduction of O3 in the stratosphere follows the descent of SPE-enhanced NO x , 

reaching losses of 500 ppbv during December and early January. Reductions of ozone 
remain at 400 to 500 ppbv (5 to 10%) from 25 to 30 km until spring. Although most of 
the ozone reduction occurs above the stratospheric ozone layer, where number densities 
peak from 15 to 20 km, ozone losses could nonetheless reduce oxidation rates in the 
425 middle to lower stratosphere and impact the chemistry, dynamics, and radiative 
properties of this region. 

Not surprisingly, the behavior of HO x , NO x , and O3 following the 9 November 2000 
SPE in these simulations is similar to the Jackman et al., [2009] WACCM results. 
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Jackman et al., [2009] show an increase in upper stratospheric NO y (mostly NO x ), 

430 exceeding 1000% for several days after the Nov 9 SPE, followed by the prolonged 

downward transport of enhanced NO y during the weeks and months following the event, 
exceeding 20% at 30 km during January and February. Model calculations of HO x and O 3 
presented in this paper are also consistent with Jackman et al. [2008; 2009]. 

Figure 5 shows a comparison of NO 2 from WACCM with POAM III (Polar Ozone 
435 and Aerosol Measurement) satellite measurements within the Arctic vortex. Descriptions 
of the POAM III data along with methods for detennining the edge of the polar vortex for 
Figure 5a are provided in Randall et al. [2002], Harvey et al. [2002], and Randall et al. 
[2007]. Scaled potential vorticity (sPV) is used to identify the location of the vortex for 
the WACCM results, consistent with the analysis throughout this paper. Magnitudes of 
440 NO 2 simulated by WACCM are similar to observations by POAM III, with background 
levels 1 to 3 ppbv between 30 and 40 km during November, increasing substantially in 
early February. The WACCM simulations do a better job capturing levels of NO 2 during 
December and January when including SPEs. Favorable comparisons for NO x and O 3 
have also been made among WACCM simulations and available satellite measurements 
445 for more recent SPEs, including MIPAS observations of the October 2003 “Halloween” 
SPEs [. Jackman et al., 2008; Jackman et al., 2009; Funke et ah, 2011] and MIPAS, MFS, 
and ACE-FTS observations for the January 2005 SPE [. Jackman et al., 2011]. 

Simulations using global models other than WACCM to study SPEs similar in 
magnitude to this study also show enhancements of NO x and reductions of O3 in the 
450 stratosphere during the months following each event [e.g., Semeniuk et al., 2005; 

Paivarinta et ah, 2013], Calisto et al. [2012] use the SOCOL (ECHAM4+MEZON) 
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model to study how a hypothetical extremely large SPE (the August 1972 SPE scaled to 
fluences associated with the 1859 Carrington Event) would impact the atmosphere if it 
occurred during a contemporary solar maximum. The short-term enhancements of HO x 
455 (10 ppbv) and NO x (200 ppbv) exceed the 9 November SPE WACCM simulations. 

However, longer-tenn increases in NO x (10%) and decreases in O 3 (20-40%) in the 
stratosphere are comparable to the 9 November SPE WACCM simulations, likely 
attributable to solar radiation variations between fall and winter. While the SOCOL 
simulation of the extreme SPE results in significant depletion of total column O3 (20 
460 DU), the WACCM simulations of the 9 November SPE show a vortex-averaged decrease 
reaching only 5 DU by late January, with reductions up to 10 DU occurring sporadically 
in the center of the vortex. 

3.3. Modeled Conversion of NO x to NO y and the Downward Transport of NO y 

465 As NO x descends within the stratosphere, oxidation reactions convert NO x into other 
NOy species, primarily through the following pathways: 

N0 2 + OH + (M) -► HN0 3 
N0 2 + N0 3 + (M) -*■ N 2 0 5 

470 (where N0 3 is produced from N0 2 reactions with O and O3) 

These products of NO x oxidation are of particular interest because of their role in nitrate 
ion deposition. Background NO y densities peak in the lower stratosphere, where a 
reservoir of HNO 3 accumulates from the oxidation of nitrous oxide (N 2 0) emitted at the 
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475 surface. The WACCM simulations presented in this paper strive to answer whether or not 
the NO y enhancements resulting from the oxidation of SPE NO x can descend and persist 
within this background pool of NO y in the lower stratosphere at levels high enough to be 
detected at the surface. Note that WACCM does not include reactions involving hydrated 
cluster ions, which have also been implicated in the production of HNO3 at altitudes over 
480 35 km [ Kawa et al., 1995; Verronen et al., 2008; Kvissel et al., 2012]. 

Figure 6 presents vortex-averaged NO y within the stratosphere from November 
through March. Figure 6a shows NO y in simulations without SPEs, depicting the 
subsidence of air within the polar vortex during winter. Figure 6b shows the 
corresponding downward transport of SPE NO y enhancements within the stratosphere. 

485 By late January there is a thin layer (~ 5 km) of 5 to 10 ppbv SPE-enhanced NO y around 
30 km. SPE enhancements continue through March, remaining at altitudes well above 20 
km. Figure 6a clearly identifies a background pool of NO y in the lower stratosphere (10 
to 15 ppbv), an order of magnitude larger than typical values in the middle stratosphere (1 
to 4 ppbv) and troposphere (less than 1 ppbv). While downward transport is best studied 
490 using mole ratios or mixing ratios (measures of composition that are independent of 
density), it is nonetheless important to draw attention to the exponential decrease in 
atmospheric density with height. A background level of 15 ppbv NO y at 50 hPa (~20 km) 
contains a factor of 10 more molecules than 15 ppbv of SPE-enhanced NO y at 5 hPa (~35 
km). Figure 6c emphasizes that NO y background number densities peak below 20 km. 

495 The challenge remains explaining how a thin SPE-produced layer averaging 5 to 10 ppbv 
NO y at 30 km might be detected at the surface given the thick background pool of 10 
to 15 ppbv NO y in the thermally stable lower stratosphere. 
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As a means of estimating the maximum potential for enhanced deposition of nitrogen 
at the surface as a result of the 9 November SPE, Figure 6d shows vortex-averaged NO y 
500 total column densities from model simulations with SPEs (blue) and without SPEs (red). 
Vertical column densities are integrated over each model surface grid, and vortex 
averages are calculated using latitudes and longitudes where sPV values on the 500 K 
isentropic surface (~20 km) are greater than 1.4x1 0‘ 4 s' 1 . Figure 6e provides similar 
model results for column densities integrated above 30 km. A large relative increase in 
505 column densities above 30 km reflects significant SPE enhancement of NO y . However, 
the absolute magnitude of the stratospheric SPE-induced NO y enhancement is small 
compared to the total column density of NO y . There is a general increase in total column 
density of NO y throughout winter. However, vortex-averaged SPE enhancements of total 
column NO y remain below 5%, with enhancements never exceeding 20% at any location 
5io within the vortex (not shown). In contrast, nitrate peaks in snow and ice show 

enhancements 4 to 5 times background levels, significantly greater than the maximum 
potential for SPE enhancement suggested by column densities of NO y . 

The fact that the increase in NO y in the atmospheric column is on average less than 
5% (and a maximum of 20% locally) challenges the theory of SPEs being responsible for 
515 spurious peaks in nitrate deposited at the surface. However, it does not rule out the 
possibility of longer-term trends in nitrate resulting from solar activity. For example, 
although Palmer et al. [2001] do not find a correlation between nitrate spikes and 
individual solar events in ice cores sampled at Law Dome, Antarctica, a statistical 
analysis reveals an 1 1% enhancement of nitrate during the 3 to 14 months following solar 
520 events, consistent with the total column nitrate enhancement calculated in this study. 
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3.4. Modeled Partitioning of NO y 

An analysis of individual NO y species within WACCM helps to identify precursors to 
nitrate ions deposited to the surface. The majority of NO y above 40 km is in the form of 
525 NO x , while HNO 3 is the principal NO y species below 30 km. Dinitrogen pentoxide 

(N 2 O 5 ) becomes significant in the 30 to 40 km layer. HO 2 NO 2 is primarily of importance 
in the troposphere. Other NO y species such as chlorine nitrate (CIONO 2 ), bromine nitrate 
(Br 0 N 02 ), and nitrate radical (NO 3 ) are small in comparison. 

Figure 7a-d presents vortex-averaged WACCM simulation results of FINO 3 including: 
530 background levels from WACCM simulations without SPEs, enhancements from SPEs, 
SPE enhancements of column densities integrated over the total atmospheric column, and 
SPE enhancements of column densities above 30 km. Throughout late December and 
most of January, SPE enhancements of FINO 3 do not exceed 1 ppbv, levels much lower 
than the 10 ppbv background pool of FINO 3 at 20 km. The largest SPE enhancements of 
535 FINO 3 occur from 20 to 30 km during late January through early March, peaking under 3 
ppbv and remaining small in comparison to background values. Vortex-averaged total 
column densities of FINO 3 show little change as a result of SPEs, and local enhancements 
only occasionally reach 2 to 3% within the polar vortex (not shown). The most significant 
SPE enhancements of FINO 3 occur above 30 km but have little impact on total column 
540 density. 

The model calculates an increase in N 2 O 5 throughout the stratosphere as a result of 
SPEs (Figures 7e-h), with most of the enhancement occurring from 30 to 40 km. In 
contrast to HNO 3 , the model shows significant enhancements (up to 40%) of vortex- 
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averaged N2O5 total column density during the months after the 9 November event. Local 
545 N2O5 total column densities are enhanced by over 100% in the center of the vortex during 
late December and January (not shown). N2O5 makes up only 5% of total column NO y . In 
order for total column N2O5 enhancements to explain nitrate spikes at the surface, there 
would need to be a mechanism involving preferential deposition of N2O5 over HNO3. To 
the authors' knowledge, no such mechanism has been suggested in the literature. In 
550 addition, Jackman et al. [2008] and Funke et al., [2011] emphasize that WACCM SPE- 
enhanced N2O5 is larger than levels from satellite measurements, hypothesizing the need 
to include water cluster ion reactions in WACCM to convert NO3 to HNO3 in order to 
simulate lower N2O5 values [e.g., Solomon et al., 1981; Lopez-Puertes et al., 2005]. This 
additional HNO3, however, would not significantly impact the variability of nitrate at the 
555 surface, since N2O5 comprises such as small percentage of total column NO y . 

3.5. Modeled Deposition of NO y 

Several mechanisms allow nitrogen in the atmosphere to deposit as nitrate ions to 
surface snow. Dry deposition involves gas phase and particulate nitrogen species sticking 
560 to the surface snow in the absence of precipitation. Wet deposition removes nitrate as 

HNO3 (and to a lesser extent HCLNCL) is taken up by falling snow by means of in-cloud 
nucleation scavenging and below-cloud impaction scavenging. Nitrate deposition can 
also result from the uptake of HN0 3 by polar stratospheric clouds (PSCs), heterogeneous 
conversion of N 2 0 5 to HN0 3 on PSCs, and subsequent gravitational settling of cloud 


565 particles to the troposphere. 
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Although background levels of NO y are large in the lower stratosphere, thermal 
stability inhibits mixing across the tropopause. In the absence of stratospheric intrusions, 
denitrification by gravitational settling of PSC particles is the most likely mechanism for 
upper atmospheric nitrate to reach the troposphere. WACCM relies on an equilibrium 
570 approach to calculate the amount of HNO3 condensed on polar stratospheric clouds, 

which in turn determines the radii and sehling velocities of these aerosols (see supporting 
information in Kinnison et al. [2007]). Although heterogeneous reactions in WACCM 
occur on PSCs composed of supercooled ternary solutions (STS), nitric acid trihiydrate 
(NAT), and water-ice [Kinnison et al, 2007; Lamarque et al., 2012; Wegner et al., 2013], 
575 condensed nitrogen is limited to STS and NAT aerosols in WACCM, with gravitational 
settling occurring only for NAT particles because of their larger radii. Once settling and 
advection bring condensed HNO3 to the troposphere, it may be removed from the 
modeled atmosphere via wet deposition. 

Denitrification through cloud sedimentation is not as common in the Arctic as in the 
580 Antarctic, where temperatures are low enough to allow PSC particles to grow to larger 
masses. Although PSCs and denitrification have been observed in the Arctic during years 
when the polar vortex is strong [e.g, Waibel et al., 1999; Kondo et al., 2000; Santee et al., 
2000; Fahey et al., 2001; Popp et al, 2001], the winter of2000-2001 exhibits neither a 
consistently strong vortex nor low temperatures (see Figure 1 from Manney et al. [2006]), 
585 In addition, observations show significant re -nitrification in the Arctic as increasing 

temperatures in the lower stratosphere and troposphere cause nitrogen from PSCs to re- 
enter the gas phase [Dibb et al., 2006]. 
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Nonetheless, for completeness, it is useful to consider model-calculated PSCs when 
examining maximum potential nitrate deposition from the stratosphere. WACCM 
590 calculates the presence of all three types of PSCs over the Arctic during winter 2000- 
2001, with STS and NAT particles most abundant in mid-January and ice clouds briefly 
present when stratospheric temperatures reach a minimum. Figure 8 presents vortex- 
averaged vertical temperatures, vortex-averaged condensed nitrate as NAT particles (the 
most relevant PSC leading to denitrification in WACCM), and condensed nitrate as NAT 
595 PSCs for 5 January 2001 (a representative day when PSCs are significant). Although 

results are given for simulations with SPEs, differences in temperature from simulations 
without SPEs are negligible and can be attributed to computational noise. WACCM 
temperatures drop below 200 K during January, reaching a minimum of 194 K at 20 km . 

Condensed NAT in the form of PSCs peaks above 15 km during January, contributing 
600 50% of total column FIN03 on 5 January 2001. However, NAT radii are small (median 

radii less than 1.5 microns) and the modeled flux of nitrate never extends below 17 km, 
as temperatures exceed 200 K and condensed nitrate ions re-enter the gas phase before 
reaching the troposphere. The impact of solar proton events on simulated vortex-averaged 
total column condensed NAT as PSCs is less than 0.2%, with local maximum 
605 enhancements of 1 to 3%. The limited impact of SPEs on nitrate in the form of NAT 

PSCs is consistent with the limited SPE enhancement of total HNO3 presented in Figure 
7. As a comparison to Antarctic simulations where temperatures are lower, NAT and ice 
PSCs are prevalent, and stratospheric denitrification has been observed to be significant, 
Howeverm Jackman et al. [1990] calculate only a 10.6% maximum SPE enhancement of 
6io nitrate deposition following the August 1972 SPE event over Antarctica. In summary, it 
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is unlikely that SPE-enhanced denitrification could account for the 4 to 5 fold nitrate 
spikes observed in Greenland surface snow. 

Enhanced nitrate from the lower stratosphere descends into the troposphere not only 
through the gravitational settling of PSCs but also from dynamics associated with 
615 stratosphere-troposphere exchange. Given the minimal SPE enhancement of NO y in the 
lower stratosphere, downward transport of SPE-enhanced NO y through isolated 
tropopause folds and stratospheric intrusions could still not explain the magnitude of 
observed nitrate ion spikes. 

In spite of WACCM’s inability to calculate enhancements of NO y from the 9 
620 November 2000 SPE large enough to explain surface spikes in nitrate ions, it is 

nonetheless worthwhile to examine how closely WACCM’s calculations of nitrate 
deposition resemble observations. Summertime estimates of the NO3' inventory presented 
by Bergin et al. [1995] attribute 93% of deposition to snow, 6% to fog and 1% to dry 
deposition. The strong stability of the wintertime surface layer would suggest that dry 
625 deposition plays an even less significant role during the time period considered in this 
study. 

WACCM treats wet deposition in the troposphere as a first order loss rate from the gas 
phase, with precipitation rates and cloud cover derived from MERRA meteorological 
fields. During each time step, WACCM calculates how much of a given gas will be 
630 incorporated into aerosols within each grid box and removes this amount from the grid 
box concentration that is passed to the next time step. While re-evaporation and 
desorption are included within a given time step, WACCM does not archive and 
propagate the soluble ion species nor account for subsequent aqueous reactions. The 
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following analysis assumes that all loss from the gas phase in the tropospheric grid 
635 column ends up in precipitated water at the end of each half-hour time step. This analysis 
estimates the potential enhancement of nitrate at the surface by wet deposition within the 
entire vertical tropospheric column, most likely an overestimate especially when 
precipitation levels are minimal. 

Figure 9 shows estimates for nitrate deposition resulting from the wet deposition of 
640 FINO 3 in WACCM directly over Summit both with and without SPEs. The method of 
calculating modeled wet deposition instantaneously from the entire atmospheric column 
is useful for determining the maximum potential deposition of nitrate but will tend to 
over-estimates nitrate deposition, as is apparent in the different scales on the left and right 
vertical axes. This variation in magnitude is confounded by a coarse grid model 
645 resolution that is unable to represent the spatial variability caused by blowing or drifting 
snow, local accumulation and ablation, or sporadic fog and rime at the surface. Instead, 
our interest lies in identifying whether or not the model can capture the relative 
variability associated with the measured nitrate ions. Indeed, the model does capture 
relative peaks in nitrate on 23 November, 14 December, and 24 January corresponding to 
650 the nitrate ion spikes of interest in this study observed in surface snow on 22-24 

November, 13 December, and 25 January. We present results for simulations with and 
without SPEs for completeness, noting that differences are negligible (never exceeding 
0.2%), once again supporting the conclusion that enhanced nitrate concentrations on 
these days are not related to SPEs. 

655 


3.6. Alternative Explanations for Nitrate Spikes 
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The WACCM simulations provide evidence that nitrate spikes not readily accounted 
for by soluble ion correlations may nonetheless be related to tropospheric sources. Figure 
10 shows the time evolution of vertical profiles of NO x and FINO 3 above Summit from 
660 November through January. Elevated NO x , with a lifetime on the order of days, indicates 
relatively fresh pollution sources. FINO 3 , a product of NO x oxidation with a lifetime on 
the order of months, is more characteristic of aged plumes. WACCM calculates enhanced 
NO x and FINO 3 in the lower troposphere during 22-24 November. High levels of NO x in 
the middle to upper troposphere occur around 14 December. Elevated levels of both NO x 
665 and FINO 3 are present in the lower to middle troposphere on 25 January. 

Figure 1 1 presents isobaric maps of NO x and NO y from WACCM, identifying polluted 
continental plumes corresponding to vertical enhancements above Summit. Figure 11a 
depicts NO x at 800 hPa (~1.5 km) in the lower troposphere above the marine boundary 
layer one day prior to each nitrate ion spike not attributable to tropospheric sources. 

670 Figure 1 lb shows NO y plots at 500 hPa (-3.5 km) on days coinciding with these nitrate 
ion spikes. Overlaid vectors indicate the direction and intensity of winds. The WACCM 
simulations show pollution from Europe reaching Summit on 22 November, with a 
polluted plume from North America also evident to the south. The middle panels in 
Figure 1 1 provide a snapshot of a polluted plume traveling from North America at low 
675 altitudes on 13 December and aloft over Greenland on 14 December. At altitudes above 5 
km (not shown in Figure 1 1), higher wind speeds steer this plume directly over Summit, 
consistent with elevated levels of NO x in Figure 10. The simulations show transport of 
NO x from Europe on 24 January, with a broad region of enhanced NO y over Greenland 
on 25 January. These model simulations of continental plumes suggest that nitrate spikes 
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680 in surface snow during these time periods are more likely the result of continental 
anthropogenic pollution than SPEs. 

4. Conclusions 

This study screens a two-year data set of daily measurements of ions in surface snow at 
685 Summit, Greenland from 2000 to 2001 for known tropospheric sources in the search for 
evidence of solar proton events in nitrate records. WACCM modeling simulations 
examine transport, chemistry, and deposition during three specific time periods when 
correlations between nitrate and other soluble ions are inconsistent with tropospheric 
sources. The model calculations confirm that solar proton events significantly 
690 impact HO x , NO x , and O 3 levels in the mesosphere and stratosphere during the weeks and 
months after the major 9 November 2000 solar proton event. However, there is never a 
time during the simulation when SPE-enhanced NO y within the atmospheric column is 
large enough to account for the observed nitrate peaks in surface snow. 

We see no convincing evidence that SPEs are related to impulsive nitrate spikes in 
695 surface snow at Summit in the winter of 2000 to 2001 but suggest that spikes not readily 
accounted for by soluble ion correlations are the result of deposition from polluted 
plumes originating in North America and Europe. This conclusion is particularly 
compelling for recent centuries when anthropogenic emissions are capable of modifying 
atmospheric composition on a global scale. 

700 The limited SPE enhancement of total column NO y (5% vortex-averaged and 20% 

local maxima) leads us to conclude that impulsive spikes of nitrate ions at the surface are 
unlikely to result from SPEs similar in magnitude to the 9 November 2000 SPE. It would 
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be worthwhile, however, to consider how large an event would be necessary to produce 
nitrate ion levels at the surface discernible from tropospheric sources as well as the 
705 likelihood of such events given the limits of solar flare energy [Aulanier et al., 2013]. 

The association of nitrate ion levels with solar activity on centennial (Gleissberg) and 
millennial timescales also remains of interest [e.g., Motizuki et al., 2009; Traversi et al., 
2012; Ogurtsov and Oinonen, 2014]. 

A promising alternative to nitrate ions in the search for proxies for historical SPEs 
7io remains the study of cosmogenic radionuclides such as Carbon-14 ( 14 C) and Beryllium- 
10 ( 10 Be) [ Steinhilber et al., 2012]. Measurements of the cosmogenic radionuclide 
Beryllium-7 ( 7 Be) are available from Summit from 1997-1998, 2000-2002, and 2003- 
present [Dibb, 2007]. Although too short-lived to serve as a historical proxy, 
measurements and model comparisons of 7 Be following recent SPEs would provide 
715 insight into vertical transport and deposition processes, paving the way for modeling 

studies involving longer-lived cosmogenic radionuclides in the context of solar particle 
events. 

Although this study could not definitively link surface observations with solar 
energetic protons impacting the upper atmosphere, WACCM results once again point to 
720 significant impacts of SPEs on the middle and upper atmospheric concentrations of HO x , 
NO x , and 0 3 , adding to the growing collection of satellite observations and modeling 
experiments that strive to clarify perturbations in the chemistry, radiation budget, and 
dynamics of the atmosphere resulting from solar variability. The quest for an alternate 
proxy to nitrate for studying solar activity through history remains compelling, 
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725 particularly with regard to protecting technological infrastructure, understanding climate, 
and validating predictive models for space weather. 
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Figures and Figure Captions 
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Figure 1 . Schematic diagram of processes involved with nitrate deposition from solar 
energetic protons. 
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Figure 2. Daily measurements of ion content (nmol/kg) in surface snow at Summit, 
Greenland from August, 2000 to August, 2002 [Dibb et al., 2007]. Concurrent peaks of 
1185 NO 3 ' and NFt 4 + , N0 3 ‘ and S0 4 2 ', N0 3 ‘ and Na + , and N0 3 ‘ and Ca 2 + are dated by color for 
source identification. Nitrate spikes dated in black represent enhancements not readily 
attributable to tropospheric sources. Average coefficient of variation (CV avg ) indicates the 
spread among three simultaneous snow samples used to calculate each daily average, 
representing the major source of error. The lower graph identifies major SPE events in 
1190 the context of daily solar proton fluxes from GOES -8 in protons/(cm 2 -day-sr) (NOAA at 
http://www.swpc.noaa.gov/). 
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1195 Figure 3. Time evolution of the vertical structure for WACCM vortex-averaged a) HO x , 
b) NO x , and c) O 3 during the weeks following the 9 Nov 2000 event (mole ratios). The 
vertical scale is represented by potential temperature ( 0 ), left, and approximate 
geopotential height (Z), right. Top: with SPEs. Middle: no SPEs. Bottom: HO x 
enhancement (HO x with SPEs - HO x no SPEs), NO x enhancement (NO x with SPEs - 
1200 NO x no SPEs), and O3 reduction (O3 no SPEs - O3 with SPEs). 
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Figure 4. Time evolution of the vortex-averaged a) enhancement of NO x and b) reduction 
O3 from November through March (mole ratios). 

1205 



N D J F M A 




Figure 5. Comparison of WACCM NO 2 with POAM III satellite observations: a) POAM 
1210 III observations, b) WACCM with SPEs, c) WACCM without SPEs, d) WACCM SPE 
enhancements (units in ppbv). WACCM results are linearly interpolated and smoothed 
for clarity during periods when satellite measurements are outside the WACCM polar 
vortex defined according to sPV > lxlO' 4 s' 1 . 
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Figure 6. WACCM vortex-averaged NO y : a) no SPEs (mole ratios), b) SPE 
enhancements (mole ratios), c) no SPEs (number density), d) column densities 
1220 (molecules cm' 2 ) integrated throughout the total atmosphere with SPEs (blue) and no 

SPEs (red), and e) column densities (molecules cm' 2 ) integrate from 30 km to the top of 
the atmosphere with SPEs (blue) and no SPEs (red). 
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1225 

Figure 7. Time evolution of WACCM vortex-averages: a) background HNO3 (no SPEs), 
b) SPE e nh ancement of HNO3 (HNO3 with SPEs - HNO3 no SPEs), c) total column 
density of FINO3, d) column density of FINO3 above 30 km, e) background N 2 0 5 (no 
SPEs), f) SPE enhancement of N2O5 (N2O5 with SPEs - N2O5 without SPEs), g) total 
1230 column density of N2O5, h) column density of N2O5 above 30 km. 
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Figure 8. WACCM vortex-averaged vertical profiles for a) temperature (K) and b) 
condensed nitrate as nitric acid trihydrate (NAT) PSC particles (mole ratios), c) Vortex- 
averaged percent of total column NAT in the fonn of PSCs with respect to total gas plus 
condensed phase FINO3. Although results are from simulations with SPEs, differences 
with respect to no SPE simulations are negligible. 
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Figure 9. Estimates of nitrate from wet deposition at Summit, Greenland during the 2000- 
2001 winter (nmol/kg). The blue and black lines along with the left axes present the 
concentration of nitrogen resulting from the loss of FINO3 in precipitation throughout the 
1245 atmospheric column during the WACCM simulations a) with SPEs and b) without SPEs. 
The red lines and right axes refer to measurements of NO3' in daily samples of surface 
snow. 
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Figure 10. WACCM time evolution of the vertical profile of a) NO x and b) HNO 3 above 
Summit, Greenland. Recall that Summit is located 3.2 km above sea level with surface 
pressures from 660-680 hPa. Peak NO x and HNO 3 are circled during time periods where 
surface snow measurements indicate nitrate ion spikes uncorrelated with other ions. 
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Figure 1 1 . Transport of polluted continental plumes simulated by WACCM prior to 
nitrate ion spikes in Summit snow on 22-24 November, 14 December, and 25 January, a) 
1260 NO x at 800 hPa (~1.5 km) on 21 November, 13 December, and 24 January. The white 
region over Greenland indicates surface elevations above 800 hPa. b) NO y at 600 hPa 
(-3.5 km) on 22 November, 14 December, and 25 January. Wind vectors are overlaid on 
both plots to indicate the travel direction of the polluted plumes. 
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